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QHlBarchan dunes — crescent-shaped dunes that form in areas of unidirectional winds and low sand availability — commonly 

'^ display an asymmetric shape, with one limb extended downwind. Several factors have been identified as potential causes 

_j_ for barchan dune asymmetry on Earth and Mars: asymmetric bimodal wind regime, topography, infiux asymmetry 

^Sj and dune collision. However, the dynamics and potential range of barchan morphologies emerging under each specific 

scenario that leads to dune asymmetry are far from being understood. In the present work, we use dune modeling in 

order to investigate the formation and evolution of asymmetric barchans. We find that a bimodal wind regime causes 

Jimb extension when the divergence angle between primary and secondary winds is larger than 90°, whereas the extended 

(y;^ limb evolves into a self dune if the ratio between secondary and primary transport rates is larger than 25%. Calculations 

^ of dune formation on an inclined surface under constant wind direction also lead to barchan asymmetry, however no 

2 self dune is obtained from surface tilting alone. Asymmetric barchans migrating along a tilted surface move laterally, 

C with transverse migration velocity proportional to the slope of the terrain. Limb elongation induced by topography can 

occur when a barchan crosses a topographic rise. Furthermore, transient asymmetric barchan shapes with extended 

limb also emerge during collisions between dunes or due to an asymmetric infiux. Our findings can be useful for making 

quantitative inference on local wind regimes or spatial heterogeneities in transport conditions of planetary dune fields 



o 

O hosting asymmetric barchans. 
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1. Introduction 

The classical symmetric shape of barchan dunes is far 
Jrom being prevalent in nature. A wide variety of barchan 
morphologies on Earth and Mars are asymmetric, wit h 
one limb extended downwind ('Xsoarl 'l984 ijBourkel . I2Q1Q[ ). 
Since the pioneering works by Bagnold (194l|), various con- 
ceptual models have been proposed in order to explain the 
existing asymmetric dune morphologies. Dune asymme- 
try has been most predominantly attributed to asymmet- 
ric bimodal winds, topography, dune collisions or asym- 
metric sand influx. Indeed, understanding dune asymme- 
try constitutes an important issue in planetary science, as 
asymmetric dunes could potentially serve as a proxy for 
local wind regimes or variations in sand supply or topog- 
raphy. However, the significance of the different causes 
of barchan asymmetry is still poorly understood, whereas 
the potential range of dune morphologies resulting in each 
case remains unknown. Although some insights into the 



dynamics of asymmetric dunes could be gained fr om field 
monitoring within a time span of a few years (jBourke, 
120101 ). assessment of dune shape evolution is a difficult 
task owing to the long timescales involved in dune pro- 
cesses. Moreoeyer, d i verse asymmetric dune mo rphologies 
( Bourkeet all . l2004l iBourke and Goudid . l2009l ) might be 
the outcome of different factors in concurrent action, which 
poses a further challenge in the investigation of the devel- 
opmental stages of asymmetric dunes through field obser- 
vations. 



Numerical modeling — which has become indispens- 
able in the investigation of aeolian and dune processes 
on Earth and oth er planetary bodies (JBourke et al.L I2OIOI : 
iKok et al.L |2012[ ) — could provide a helpful tool in the 
study of barchan dune asymmetry. In order to shed light 
on the factors competing for the appearance of various 
asymmetric dune shapes, a minimal model that accounts 
for a mathematical description of the wind over dunes, 
as well as for the equations for grain transport and land- 
scape evolution is required. Such a mo del has been devel- 
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oped i n the cours e of the last decade (ISauermann et al 
200 ll : iKrov et al.L I2OO2I : iDuran and Herrmannl . l2006a : 
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Duran et al.L |2Q1Q| ). This model combines an analytical 
model for the average turbulent flow over the sand land- 
scape with a continuum model for sand transport. The 
model has proven to reproduce the shape of different types 
of dunes with good quantitative agreement with mea- 
surements ( Sauermann et alJj 2QQ3: D uran and HerrmannI , 
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Parteli and Herrmann . .20071 ). The model has 



12006, 

been further extended in order to model longitudinal self 
dunes and diverse un usual Martian dune forms under bi- 
modal wind regimes (jParteli et al.L l2009l ). 

In the present work, we use the dune model to investi- 
gate the role of different asymmetry causes for the diversity 
of asymmetric barchan shapes reported previou sly from 
obser vations of dune fields on Earth and Mars (jBourkd . 
20101 ) . We calculate the evolution of a symmetric barchan 
dune under the separate action of the potential causes for 
dune asymmetry. 

This paper is organized as follows. In Section [2] the 
dune model is described. In Sections [3l|6] we present and 
discuss the results obtained from calculations of barchan 
dune asymmetry due to bimodal wind regimes, topogra- 
phy, asymmetric sand supply and dune collision, respec- 
tively. Conclusions are presented in Section [71 



2. Model description 

The model used in the calculations of the present 
work consists of a set of mathematical equations that 
compute: the average turbulent wind field (the sur- 
face shear stress (r) over the topography; the mass flux 
(q) of saltating particles due to the shear stress, and 
the time evol ution of the surface resulting from parti- 
cle transport ([Sauerniann et al.L uOOlc iKrov et al.L l2002l : 
Duran and Herrmannl . l2006al ). In this Section, we present 
a brief description of the model equations and of the cal- 
culation procedure. 



2.1. Wind field 

The first step consists of describing quantitatively the 
average three-dimensional turbulent wind field over the 
dune. The average shear stress field (r) is calculated 
through solving a set of analytical equations developed by 
Wenget al.l (|l99lh . 

In the absence of dunes, the wind velocity v{z) within 
the atmospheric boundary layer increases logarithmically 
with the height {z) above the flat ground, i.e. v{z) = 
2.5n*oln(z/zo), where the wind shear velocity u^o is 
used to define the (undisturbed) shear stress tq = 
Pfluid|'^*o|'^*05 and p^md is the air density. A smooth 
hill or dune introduces a perturbation in the wind field. 
The Fourier-transformed longitudinal and transverse com- 
ponents of the shear stress perturbation (f) due to the 
local topography are computed using following equations. 
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where x and y are parallel, respectively, perpendicular to 
the wind direction, a = y^iLkxZo/l^ kx and ky are the 
components of the wave vector ^, i.e. the coordinates 
in Fourier space, Kq and Ki are modified Bessel func- 
tions and h^ is the Fourier transform of the height profile; 
U is the vertical velocity profile which is suitably non- 
dimensionalized, / is the depth of the inner layer of the flow 
and L is a typical length scale of the hill or dune and is 
given by 1/4 the mean wavelength of the Fourier represen- 
tation of the height profile. In the presence of the saltating 
grains, the roughness of the dune's surface incre ases to an 
apparent value Zq, the aerodynamic roughness (^Bagnold], 
1941). The wind velocity over the dune is calculated us- 
ing Z q = 1mm, a value based on experimental observations 
(Bagnoldl. Il94ll : lAndreotti [20q3 ) . The surface shear stress 
field is obtained, then, with the equation. 



T = |to|(to/|to| +f). 



(3) 



where tq is the undisturbed shear stress over the flat 
ground. 

Flow separation at the dune brink due to the strong lo- 
cal curvature of the surface gives r ise to a zone of recircu- 
lating flow (.Walker and Nickling. . i2002l : [Herrmann et al. 



200 5), which cannot be described by the analytical model 
(|Weng et al.l Il99l[ ). This problem is overcome introduc- 
ing the so-called "separation bubble" : for each longitudi- 
nal slice of the dune, a separation streamline connecting 
the brink to the reattachment point is introduced at the 
lee. Each separation streamline is fitted to a third-order 
polynomial, the par ameters of which ar e determined as de- 
scribed in detail by iKrov et al.l (|2002[ ). The wind model 
( Weng et al.l Il991[ ) is then solved for the smooth "enve- 
lope" th at comprises the separation bubble and the dune 
surface rKrov et al.Ll2002[ ). Thereafter, the shear stress in 
the recirculating zone within the separation bubble is set 
as zero, since the net transport within the bubble essen- 
tially vanishes. 

2.2. Sand flux 

Next, the mass flux of particles in saltation — which 
is the dominant transport mode and consists of grains 
hoping in ballistic trajectories and ejecting new parti- 
cles upon collision with the ground — is computed. The 
saltation cloud is regarded as a thin, fluid-like layer 
that can exchange san d with the immobile sand bed 



( Sauermann et al.Ll200l[ ). 



When the wind shear stress exceeds a minimal thresh- 
old and saltation begins, the sand flux first grows expo- 
nentially due to the multiplicative process inherent to the 
saltation process. However, since saltating grains acceler- 
ate at cost of aeolian momentum, the flux cannot increase 
beyond a maximal value, the so-cahed saturated sand flux 
( Bagnoldl . Il94ll : I Sauermann et al.l I2OOII : lAlmeida et al 



2QQ8I : iDuran et all l2Qlll : IPahtz et all 120121 ) which is 
reached after a saturation transient where the air shear 
stress within the saltation cloud equals the minimal value 
for susta ined saltation, i.e. the impact threshold, rt = 
Pfluid'^*t teagnoldl , Il941f ). defined in terms of the impact 
threshold shear velocity ii*f In the continuum model, the 
derivation of the three-dimensional equations for the sand 
fiux is performed by explicitly accounting for the decrease 
in aeolian shear stress due to the growth of the number of 
particles in saltation (the "feedback effect" (jOwenlJ 19641 )), 
which leads to the saturation transient of the fiux as de- 
scribed above. The following equation is obtained for the 
sand fiux {q) over the terrain. 



V-q=(l-|q|M)|q|/4 



(4) 



where ^s = [2Q^|'i^s|/5'](T — rt) is the saturated fiux; 4 = 
[2a I Vs P /gj] Tt (r — Tt ) ~ is the characteristic length of fiux 
saturation; g is gravity; a ^ 0.4 and 7 ^ 0.2 are em- 
pirica l ly determined model pa r ameter s ([Sauermann et al.L 
20011 iDuran and Herrmannl . l2006a[ ). and the steady- 
state velocity of the particles in saltation (vs) is cal- 
culated numerically from the balance between aeolian 
drag, gravitational and bed friction forces on the parti- 
cles (|Duran et al.L I2OIOI ) . The outcome of the calculation 
is a two-dimensional field {q{x^y)) that gives the height- 
integrated, average mass fiux of saltating particles over the 
terrain. 



2.3. Surface evolution 

Finally, mass conservation is used in order to compute 
the evolution of the local height (/i(x, y)) through the equa- 
tion, 

Ps^nddh/dt = -V • gr, (5) 

where Psand is the bulk sand density. Eq. (jFj) implies that 
the deposition (erosion) occurs at those places where the 
fiux locally decreases (increases) downwind. 

Avalanches — wherever the local slope exceeds the an- 
gle of repose of the sand {0^ ^ 34°), the surface is re- 
laxed through avalanches in the direction of the steepest 
descent. Avalanches are considered to be instantaneous 
as their time-scale can be regarded as negligible compared 
to the time-scale of the surface evolution due to aeolian 
transport. The fiux of avalanches along the slip- face is 
given by. 



<7avai = ^[tanh(V/i) - tanh(6>dyn) 



Vh 



(6) 



where k = 0.9 is a parameter and ^dyn = 33° is the so- 
called "dynamic" angle of repose, wh i ch cha racterizes the 
surface after relaxation ( Duran et al.L I2OIOI ). The update 
in the local height is obtained by solving Eq. (|5|) using the 
fiux due to avalanches given by Eq. ([6j). The calculation 
is repeated until the local slope is below ^dyn- 



2.4- Description of the calculation procedure 

Summarizing, the model consists of iteratively perform- 
ing the following calculations: 

1. the average shear stress (r) over the surface is com- 
puted using Eqs. ([I]), (|2|) and ^] 

2. next, the height-integrated average mass fiux q over 
the terrain is calculated by solving Eq. (JH); 

3. the change in the local surface is computed with Eq. 
dU; wherever the local inclination is larger than 34°, 
the fiux due to avalanches is calculated using Eq. (J6j) 
and the topography updated again using Eq. (|5|). 

The initial surface is a smooth hill of Gaussian shape, 
which is subjected to a constant upwind shear stress of 
value To — in the following, we call the upwind shear stress 
simply r. Calculations are performed with open bound- 
aries and an infiux q^^ at the inlet, which is typ ically 20% 
of the saturated fiux, q^ (jFrvberger et al.L Il984l ). 



3. Bimodal wind regimes 

Most attempts to model barchan asymmetry have con- 
centrated on the role of asymmetric bimodal wind regimes 
for the elongation of one barchan limb. The bimodal wind 
is said to be asymmetric when the transport rates of both 
wind components are not equal. 

The first conceptual model was by Bagnoldl (|194l l). Ac- 
cording to this model, a symmetric barchan, originally 
formed by a gentle wind, becomes asymmetric if a storm 
wind blows from a secondary direction, making an ac- 
cute divergence angle with the primary one. The limb 
exposed to the storm wind elongates as it enters the sand 
stream of the limb at the opposite side, th us evolving 



into a longi t udinal self dune (McKee and TibbittsI, 
T^QarLll982Lll983LlBristow et al.Ll2000LlRubin eraL 



1964 : 



2008; 



Tokand . |2010[ ). According to Bagnold's model, the self 
dune formed in this manner al igns approxima tely parallel 
to the direction of the storm ( Bagnoldl , Il94ll ). This con- 
ceptual mod e l was referred to by n iany authors in the past 
(IVerstappen . 1968L Ruhe . 1975 : Goudie and Wilkinson 



I977I : iMabbuttl . ll977L iLancasterL Il980[ ). though reports o n 
supporting field examples were scarce (Lancaster, 1980|). 
A different model was proposed later by Tsoar (1984): a 
gentle wind (not a storm wind) blows from the secondary 
direction, and the limb that elongates is the one opposite 
to the secondary gentle wind (Tsoar, 1984). Some field o b- 
servations exist that support this model ([BourkeL |2010[ ). 

In our calculations, we simulate a bimodal wind by 
periodically alternating the orientation of the field be- 
tween two d i rectio ns forming a divergence angle, 0^ 
( Par teli et al.L |2009[ ) — such as in the turntable experi- 
ments of ripples and subaqueous dune formation on a sed- 
imen t bed under bimodal fiow regimes (|Rubin and Hunterl , 
19871 : iRubinand Ikeda|, Il990l ). The wind model is solved 
considering a constant wind blowing over the rotated land- 
scape, while the separation bubble adapts t o the wind di- 
rection following the rotation of the field (jParteli et al 




Figure 1: Elongation of one barchan limb due to a bi- 
modal wind regime with obtuse divergence angle. Wind 
directions are indicated by the arrows (the primary wind, 
which has duration T^i, blows from the left). The barchan 
shape depends on r = Tw2/^wi, i-e. the duration Tw2 of 
the secondary wind relative to T^i- In these simulations, 
Twi is about 3% of the migration time of the barchan dune, 
Tm- The shapes shown for r < 25% are stationary shapes, 
whereas for r > 25%, the elongated dune horn increases 
with time. Simulation snapshots correspond to time t of 
the order of 5Tm. The barchan obtained in the simula- 
tion with r = 5% has length and width of approximately 
120 m. 



20091 ). The primary wind direction, i.e. the one that forms 
the barchan, has duration T^i and upwind shear stress ri. 
The secondary wind direction makes an angle 0^ with the 
primary wind, and lasts for a time Tw2 with upwind shear 
stress r2. 

We begin our study with an obtuse divergence angle 
(6>w = 120°), and consider two models: 

Model #1: both wind directions have the same upwind 
shear stress (ri = T2), however the primary wind has a 
longer duration (T^i > ^^2) — sl scen a rio r eminiscent 
of the experiments by iRubin and Hunter! (|1987[ ) with the 
difference, however, that in the simulations the ground is 
not covered with sand. Figure [T] shows snapshots of the 
dune evolution obtained from the calculations. As can 
be seen from Fig. (TJ the barchan develops an asymmetric 
shape, which depends on the ratio r = Tw2/^wi- If r > 
25%, then the limb at the side opposite to the secondary 
wind elongates into the resultant wind direction (Fig. [1]) 
to form a longitudinal dune. It is interesting that the 
longitudinal alignment also appears to be lost in turntable 
experiments of ripples on a sand bed when the transport 
rates differ by a factor larger than 4 (JRubin and Hunter . 
19871 ). When r = 1.0, the ba rchan dune g i ves pl ace to a 



symmetric longitudinal dune ([Parteli et al.L l2009l ) . 

Model #2: the primary wind has a la rger shear str ess 
(ri > T2) — as in the conceptual model bv lTsoarl (|l984l ) — 
while Twi = Tw2. In this case, the dune shapes do not dif- 
fer much from the ones in Fig. [H The limb opposite to the 
secondary wind elongates when the relative values of bulk 
sand flux in the secondary and primary wind directions. 



respectively Q2 and Qi, are such that r = Q2/Q1 > 25%. 
Indeed, we also considered the situation where the storm 
wi nd is the seconda ry one {Q2 > Qi), such as in the model 
bv iBagnoldl (|194ll ). In this case, primary and secondary 



winds simply exchange their roles, whereas the same dune 
shapes such as in Fig. [T] are obtained. In summary, Bag- 
nold's hypothesis that the elongating limb is the one ex- 
posed to a secondary storm wind is not supported by the 
simulation results. 

In a general manner, the condition for elongation of the 
asymmetric limb due to a bimodal wind of obtuse diver- 
gence angle, as found from our calculations, reads: 



[Q2-Tw2]/[Qi-Twi]>l/4, 



(7) 



where the indices 1 and 2 refer to the primary and to the 
secondary wind, respectively. Equation ([7]) allows us to es- 
timate the onset for emergence of asymmetry in barchan 
dunes, where both the relative duration and shear stress 
values of the wind components are accounted for quanti- 
tatively. 

Next, we extend the calculations to different values of 
the divergence angle 6^. In Fig. [2l the dune shape emerg- 
ing from an asymmetric bimodal wind with r = 50% is 
shown as a function of 0^ and ti = T^i/Tm, which is 
the duration of the primary wind, T^i, rescaled by the 
migration or reconstitution tim e of the barchan, T^ ^ 
0.02WVQi (JDuranet al.Ll2010[ ). The dune turnover time 
is roughly the time neede d for t he dune to cover a dis 
tance of its width (W) (lAllenl. 1 19741: [Lancaster . 
Rubin and Ikedal . Il990l : iHersen et al.L l2004l ). 



19881 : 



As can be seen from Fig. [2l dunes orient transversely or 
longitudinally to the resultant transport trend, depending 
on whether 6^ is acute or obtuse, respectively. Bimodal 
wind regimes with 6^ < 90° lead to rounded (or oblate) 
barchans, whereas larger values of ^w yield asymmetric 
dunes tha t elongate in the resultant transport direction 
( Parteli et al. . 2009) . However, the elongating limb is not 
stable if 0^ is within the range 90° < 6>w < 112°. In 
this case, if the influx is sufficiently small, the extended 
limb se parates from the barchan , giving rise to a ^ ^ mixed - 
state" (IRirbin and Huiiterl . Il987l : IHubin and Ikedal . Il990l : 
Parteli et al.L ^2009*), as depicted in Fig. [3^1. In fact, 
mixed states analogous to the ones found for asymmetric 
barchans have been also found in turntable experiments 
of aeolian ripples and subaque ous bedforms on a sand 
bed for the same rang e of 0^ (JRubin and Hunterl . 119871 : 
Rubin and Ikedal . Il990[ ). At very large divergence angles, 
the resulting mor phology resemb l es an asymmetric dune 
of "slim" shape ( Long and Shard . Il964t iBo urke. 2010), as 
shown in the diagram of Fig. [2] and also illustrated with 
the calculation using 0^^ = 170° in Fig. |3t. The dunes in 
Fig. [3t are "reversing" dunes, which are formed by winds 
coming from opposing dire ctions. This dune form has been 
discussed by Tsoarl (|200ll ). who also gave examples of re- 
versing dunes occuring in South Africa (Tsoar, 2001). 

Asymmetric barchans such as the ones in Fig. [1] form 
when the relative duration of the primary wind {ti = 



10° C ^ 
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Figure 2: Dune shape as a function of the divergence an- 
gle (6>w) of the bimodal wind, and of the rescaled duration 
of the primary wind (t\ = Twi/Tm), for r = 50%. The 
directions of primary and secondary winds are indicated 
by the arrows (primary wind blows to the right). Sim- 
ulation snapshots correspond to time t ^ 3Tm, except 
those obtained with ti = 10~^, which correspond to time 
t ^ 0.2T^. 
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Figure 3: Time evolution of the barchan shape under 
bimodal wind regimes with obtuse divergence angle and 
r = 50% (c.f. Fig. [2]). Time is in units of the turnover 
time of the barchan (Tm). (a) The divergence angle is 
^w = 100°, the duration of the primary wind (to the right) 
relative to T^ is ti =0.1, and the influx (gin) is 10% (top) 
and 20% (bottom); (b) l9w = 120°, gin = 20%, ix = 1.0 
(top) and 0.1 (bottom); (c) 0^ = 170°, ^in = 20%, h = 0.1. 



Twi/Tm) is within the range 10"^ < h < 10"^ When 
ti is smaller than about 0.1%, no limb elongation occurs; 
both limbs, as well as the slip face, disappear and a dome- 
like shape is obtained (c.f. Fig.[2j). As ti approaches unity, 
any change due to the secondary wind is fully compensated 
by the primary wind during every cycle of the bimodal 
wind regime — the resulting morphology is essentially a 
barchan that is periodically reformed by a unimodal wind 
of strength ri, as can be seen from Fig. [2j The dune 
shapes s hown in Fig, p for ti = 1.0 are thus "transitional" 
shapes (|de Horn l2006h . The oblique incidence of the sec- 
ondary wind leads to destabilization of the dune surface 
and the emergence of smaller barchans, which detach from 
the larger dune and migrate on the bedrock thereby alter- 
nating between both directions of the bimodal wind. Fig- 
ure [3)3 shows the calculation of a barchan under asymmet- 
ric bimodal wind regime with 6^ = 120° for ti = 1.0 (top) 
and ti =0.1 (bottom). On the basis of Fig. [SJd, it is pos- 
sible to understand the coexistence of barchans and other 
complex bedforms shaped by multimodal wind regimes. 
If the dune is small enough, it has a small turnover time 
(~ large ti) and can readapt to the prevaihng transport 
trend notw i thsta nding the complexity of the wind system 
(JLancastedJlQSSh . 



4. Topography 

Field evidence of b archan limb elongation caused by to- 
pography are uncle ar ( Finkeill959l : lLong and SharplJl964 : 
ICavl , I1999I : iBourkd , |201Q| ) . We investigate two cases where 
topography is observed to trigger dune asymmetry under 
both constant wind direction and sediment supply. 

Sloping terrain — We consider a barchan dune migrat- 
ing on a sloping terrain that makes an angle (f with the 
horizontal (Fig.|4^). The simulation starts with a symmet- 
ric barchan of height H and width W^ which was shaped 
on a flat surface under constant upwind shear stress r and 
constant influx {qin/Qs = 0.2). Then, the barchan is let 
to evolve on the inclined surface. Due to gravity, sand 
transport on the tilted surface has now a component per- 
pendicular to wind trend. 

Figure [3)3 shows that the emerging barchan shape is 
asymmetric, yet none of its limbs elongate to form a longi- 
tudinal dune dParteh et al.L 120091 : iReffet et al.Ll2010l ). The 
asymmetry arises from the combined effect of both down- 
wind sand transport and gravity-driven mass flow in the 
direction orthogonal to the wind. Sand transport along 
the lower limb has a small net component downhill due 
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Figure 4: (a) Barchan dune of cross-wind width W ^ 90 m migrating along a sloping terrain. The surface is tilted by 
an angle if. Arrows indicate the wind direction (r is the upwind shear stress) and gravity {g). (b) Dune shape as a 
function of if and r/rt (wind blows from left to right). Simulations were performed with gin/^s = 0.2. (c) Main plot: 
Ratio between transverse and longitudinal migration velocities, respectively ^t and ^l- The best fit to the data using 
the equation ^t/'^l = kqtdJiip gives kq ^ 0.26 (continuous line). Inset: kq as a function of Qin/Qs- (d) Main plot: The 
degree of asymmetry is quantified in terms of the ratio w/W ^ calculated as a function of the surface inclination {(p) and 
for the same values of r/rt shown in the main plot of (a). Inset: Length-to- width ratio {L/W) as a function of r/rt 
calculated for different values of ip within the range < tan(/? < 0.4 and for two values of (^ above this range ((/p = 0.45 
and ip = 0.50). The best fit to the data within the range < id^nip < 0.4 using the equation L/W ~ A[r/rt — 1]~ -\- B 
gives A ^ 0.42 and B = 0.8 (continuous line). 



to the tilting, and so this limb stretches downwhill. In 
contrast, the upper limb is not strechted because the sand 
transported laterally due to the tilting is trapped at the 
slip face. 

A consequence of the surface inclination is the downhill 
extension of the barchan. The migration velocity of the 
barchan dune orthogonal to the wind trend increases with 
(/?, as shown in the main plot of fig. Ht. For moderate 



values of cp within the range < tdiUip < 0.4, the ratio 
between transverse and longitudinal migration velocities, 
^T and vl, respectively, can be described by the equation: 



Vt/vi, = kqidLiiip, 



(8) 



where the coefficient kq ^ 0.26 has a slight dependence on 
the upwind sand ffux, Qin/Qs (cf. inset of fig. Ht). The 
ratio vt/vi, for barchans migrating on a sloping terrain is 



thus a function of the slope (/?, being independent of the 
wind velocity. 

The topography-induced barchan asymmetry can be 
quantified in terms of the relative width w/W^ where W 
is the total dune width, and the width w is the distance 
— measured in the direction orthogonal to wind trend — 
between the border of the dune at the lowest elevation 
and the dune's windward foot (c.f. fig. |3]i). If there is 
no surface tilting ((/? = 0), then w = 0.5 since the dune is 
symmetric along its central axis. As cp increases within the 
range < tdiUif < 0.4, w increases roughly linearly with 
tan(/?, as shown in fig. [Hi. Within this moderate range of 
(/?, the ratio w/W can be described by the equation. 



w/W = kujtdinip, 



(9) 



where the value k^ ~ 0.33 is obtained from the best fit 
to the simulation data (c.f. Fig. |3]i). The aspect ratio 
of the dune (L/W \ which is a function of the relativ e 
shear stress (r/rt) (JKrov et all l2002l : IParteH et al.LlioQTl ). 
is essentially independent of cp. In other words, r/rt and 
(f can be approximately estimated from L/W and w/W 
of the asymmetric barchan shape, respectively, provided 
other asymmetry causes (i.e. bimodal wind, dune colli- 
sions or infiux asymmetry) are not relevant at the dune 
field considered. 

Topographic rise — Asymmetric barchans occasionally 
occur approac h ing a nd crossing a topographic break in 
slope (JBourkd , 120101 ). We investigate, using the dune 
model, the evolution of a barchan crossing a ridge that is 
placed obliquely to the direction of motion (Fig. [5j). The 
dune is subjected to an unimodal wind regime and an in- 
fiux that is 20% of the saturated fiux. We consider that 
the longitudinal axis of the ridge forms a moderate an- 
gle with the direction orthogonal to the wind (about 45°). 
The ridge has Gaussian cross section, width Lr and height 
Hr, while the barchan has a height H ^ 9 m. 

Indeed, the height-to- width ratio Hr/Lr (or aspect ra- 
tio) of the ridge cannot be too large, as the ridge would act 
as a barrier for sand transport. On the contrary, a ridge 
of too small aspect ratio has negligible effect on the dune 
shape. An asymmetric dune shape with limb elongation 
occurs when the aspect ratio of the ridge is about 1/20. 
In Figs. [5)3 and [5t, the ridge has height Hr = 7 m and 
width Lr ~ 150 m. As the barchan approaches the ridge, 
deposition occurs at the lee of the ridge, on the side cor- 
responding to the closest limb. This limb then elongates 
once it arrives at the lee of the obstacle. The taller the 
ridge, the more accentuated the extension. As the other 
limb approaches sufficiently close to the ridge, a similar 
process occurs and thus that limb also elongates. 

The deposition at the ridge's lee occurs due to a decrease 
in wind strength just downwind of the crest. When the 
sand released through the limbs crosses the ridge's crest, 
it cannot escape the lee since the transport rate there is 
reduced. As the dune advances and more limb sand is 
deposited, the wind strength at the lee suffices to extend 
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Figure 5: The potential inffuence of topography in the 
development of barchan asymmetry in Peru, (a) Barchan 
dunefield migrating across a landscape where resistant lay- 
ers in underlying bedrock protrude as low ridges. Illumi- 
nation is from the top and north is to the top of the image. 
The box indicates a barchan that extends its eastern limb 
across the topographic rise. Wind regime in the area is bi- 
modal resulting in the extension of western limbs, except 
where infiuenced by topography. Images used with permis- 
sion from Google Earth, (b) Asymmetric barchan dune of 
height '^ 9 m produced with the model using a topographic 
rise of Gaussian cross section which has along- wind width 
200 m and height 6 m. (c) Simulation snapshots of the 
evolution of the barchan dune. Time increases from left 
to right and wind blows from the bottom. The simulation 
was performed with r/rt ^3.9 and Qin/Qs = 0.2. 



the Hmb, yet accumulation ensues. Future studies should 
focus on the asymmetric dune resulting from simulations 
using different angles between ridge orientation and wind 
direction, as well as different shapes of the ridge. Indeed, 
the results of our calculations should be valid for ridges 
which have a smooth cross section. If the ridge's cross 
section has a sharp edge, then ffow separation occurs and 
the ffow is diverted parallel to the ridge along the lee side 
( Tsoarl . [2001). Due to the assumption of our model that 
net fiow ceases within the separation zone, the sand de- 
posited at the lee of a ridge with sharp edge cannot be 
transported at all. Thus, three-dimensional flow patterns 
at the lee of dunes or obstacles should be accounted for 
in the modeling of dune asymmetry when the topographic 
rise has a sharp slope (Tsoar , 200 1). 

The asymmetric barchans of the field example in the 
image in Fig. [5^ display superimposed substructures or 
"sand-wave instabilities" ( Elbelrhiti et al. , 2005) indicat- 
ing the presence of a secondary wind direction. Certainly, 
the presence of oblique winds play an important role for 
the shape of the elongating limb and of the barchan. How- 
ever, our calculations show that such secondary winds are 
not a necessary ingredient for limb elongation, i.e. the 
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Figure 6: (a) Time evolution of a barchan of width 80 m subject to an asymmetric influx: The influx is qi = 20% of q^ 
in the white areas, in the dark areas the influx is ^2, where ^2/^5 = 90%, 80% and 70%, respectively from top to bottom. 
The wind blows from the left with upwind shear stress r = 3.24rt. Snapshots at time interval At = 3Tm where T^ is 
the turnover time of the barchan at t = 0. (b) k; is the partial fraction of the barchan that is within the stream of larger 
influx ((72 = 0.7, t = 22 Tm), whereas W is the total width of the barchan. (c) Main plot: w/W as a function of time for 
Q2/Qs = 0.4 (squares), 0.6 (triangles) and 0.9 (circles). The straight lines are the best fits to the simulation data using 
a linear relation w/W = kt/T^, where k = vx/^m and vt is the migration velocity in the direction orthogonal to wind. 
The inset shows that '^t is less than 10% of '^m- 



asymmetric dune pattern in Fig. [5] can be explained only 
on the basis of a variation in wind strength induced by 
a topographic obstacle placed obliquely to the direction 
of motion. The asymmetric shape of dunes near a topo- 
graphic rise can also be modified due to the presence of 
other dunes nearby, since different types of barchan asym- 
metry are generated by asymmetric influx or collisions be- 
tween dunes, as will be discussed next. 

5. Asymmetric sand supply 

The role of asymmetric sediment supply for barchan 
asymmetry has re mained uncertain in th e few i nvestiga- 
tions undertaken (|Rim|, Il958l : iLancasterl , 1 19821 : iBourkd . 
2010). It is admittedly a difficult task to conduct system- 



atic field studies on the barchan shape as a function of the 
degree of spatial asymmetry in the upwind sediment bud- 
get. The sand flux onto the windward side of a barchan is 
largely dependent upon the spatial distribution of the up- 
win d dunes or other sand sources, e.g. at a sand sea mar- 
gin ([Lancasterl , Il982[ ) . Spatial inhomogeneities in the local 
physical properties of the interdune terr ain may also play 
an im portant role for the incoming flux (|Frvberger et al. , 
Il984[ ). Furthermore, the shape of a barchan migrating in 



a fiel d is inevitably a f fected by wind trend variations in 
time (JElbelrhiti et al.l |2005[ ) , as well as by the interaction 
with neighbouring dunes through merging, breeding and 
lateral linking. 

With the help of simulations, the barchan shape re- 
sulting from asymmetric upwind flux alone can be stud- 
ied without regard to other physical factors. We con- 
sider a symmetric barchan initially migrating within an 
area where the influx qi is approximately a constant frac- 
tion (^ 20%) of the saturated flux (q^) — the value 
of qi/qs ^ 0.2 is, in fact, within average values mea- 
sure d at interdu ne areas in typical barchan dune fields 
(Frvberger et al.L 1984). The barchan has height H ^ 6.0 
m and width 80 m, and has been shaped by a constant 
shear stress r ^ 3.24rt. 

At t = 0, the upwind influx becomes asymmetric. The 
influx on the side of the barchan's left limb has a larger 
value, q2 > qi (dark areas in Fig. [6^). This model for 
asymmetric influx constitutes a simple model for a sce- 
nario typical for real dune fields, where in some areas the 
presence of a dome, a sand sheet or a larger dune migrating 
upwind may significantly increase the local interdune flux. 
Figure [6^ shows different snapshots of the evolution of the 
barchan dune subject to the asymmetric influx. Because 



more sand deposits in the areas of larger influx (^2), the 
dune adapts to the asymmetric influx by moving sideways 
towards the area where the influx is q2. Different (tran- 
sient) asymmetric dune shapes appear as the dune grows 
and moves sideways, as can be seen in Fig. [6^. The shorter 
limb is the one exposed to the larger influx value, because 
a large upwind flux prevents erosion at the windward foot 
and downwind motion of the dune. Thus, the barchan 
limb on the side of l ower influx appears advanced in rela- 
tion to its — "fat" (jLong and SharpI , Il964l : IParteli et al. , 
20071 ) — counterpart. 



most entirely into the region of larger influx and is nearly 
symmetric again. This situation is indeed observed in the 
simulation (see e.g. Fig. [6^ for q2/(ls = 0.9). However, 
as the influx becomes very large the barchan shape gives 
place to a flat dome-like dune without slip- fa ce, which may 



serve as source of sand for smaller dunes ( Duran et al. 



The fraction of the dune width {w/W ^ c.f. Fig. [Gb) that 
is within the area of larger influx ((72) increases in time, as 
can be seen in Fig.[6t. The rate with which w/W grows in 
time is deflned as the lateral migration velocity {v^) of the 
dune due to the asymmetric influx. Firstly, w/W increases 
roughly as a linear function of time which means that in 
this initial stage of lateral migration ^'t is approximately 
constant for given values of gi and (72 relative to q^. The in- 
set of Fig.[6t shows that v^ is indeed much smaller than the 
longitudinal migration velocity of the barchan (vm) even 
for large ^2/^5- Indeed, the main plot of Fig. [6t suggests 
that, after a sufficiently long time, w/W asymptotically 
approaches a maximal value (and thus v^ -^ 0). Since 
the increase of w/W with time becomes extremely slow 
after long time, and since the dune is increasing rapidly 
in volume as ^2 approaches q^^ it is difficult to reliably 
estimate the steady-state value of w/W through numeri- 
cal simulations. However, we can approximately estimate 
the steady-state value of w/W (which we call Woo/W by 
means of a simple calculation, as we will explain in the 
next paragraph. 

A barchan under constant influx q^j^ is an unstable ob- 
ject. While the dune is gaining sand due to g'ln, it is also 
losing mass through the limbs. The mass balance deter- 
mines whether the dune grows or shr inks with time. As 
shown previously (|Duran et al.L [2010), there is a critical 



2OIOI : iLuna et all , l201li I2OI2I ). We note that an increase 
of gi implies a decrease in the influx asymmetry, since q2 
is within the range qi < q2 < <7s- Thus, as qi increases, 
the steady-state shape of the dune becomes less symmetric 
(though the barchan shape gives place to a dome as dis- 
cussed above) and Woo/W ^0.5. Evidently, for qi = q2 
there is no asymmetry in the influx and the dune shape is 
thus symmetric {woo is equal to W/2). 
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Figure 7: Steady-state value Woo/W of the fraction w/W 
of the dune cross-wind width within the region of larger 
influx (^2), as a function of q2/qs for different values of 

Qi/Qs' 



value of (7in, denoted by q^, above which the dune volume ^ Dune collisions 
increases in time and below which the dune shrinks. This 
critical influx is approximately equ al to 18% of the sat- 



Dune collisions are well-known factor for the emer- 



urated flux, q^ ([Puran et all |2010[ ). The rate dW/dt at gence of comp l ex, as ymmetric dune morphologies (Tsoar, 



which the cross-wind width W of the barch an changes in 
time approximately scales with [q^^ — qc]/W (|Duran et al.L 
20101). Using this relation, we can approximately estimate 
the growth rate of w and oi W — w^ which denote the 
fractions of the dune width under influx values q2 and qi , 
respectively (cf. Fig. [GJd). The growth rates dw/dt and 
d{W — w)/dt should be proportional to [^2 — ^0]/'^ and 
[^1 — (lc]/iW — w)^ respectively. The steady-state value of 
w^ i.e. Woo is achieved when both rates equal, which gives. 



w 



Q2 



(10) 



qi^q2-'2-qc 

with gc ^ O.lSg's as mentioned above. Figure [3 shows how 
'^00/1^ depends on ^2/^5 for different values of gi/^s, with 
Qi^Q2 > Qc- For values of qi close to ^c, Woo/W approaches 
unity as (72 increases, which means that the dune enters al- 



2OOII : iBourkd l2010f ). The dynamics of collisions be- 
tween two barchans have been systematically inves- 
tigated in previous studie s using the present model 
(Schwam mle and Herrmanm 120031 : iDuran et all l2005l 
2009). It was shown that, if the symmetry axes of both 
dunes are aligned, then the collision dynamics is deter- 
mined by the volume ratio of the dunes, whereas the result- 
ing morphology is symmetric (Duran et al. , 2005) . How- 
ever, if the collision occurs with a later al offset (cf. Fig. [8]) , 
asymmetric dune shapes can arise (iCl ose-Arc edua 1 19691 : 
Hersen and Douadvl , 120051 : iDuran et al. . 2009) . 

Figure [9] shows snapshots of calculations using differ- 
ent values of volume ratio and lateral offset — deflned as 
A = \Y — y\/{W/2)^ where Y (y) is the crest's position 
of the large (small) dune in the direction orthogonal to 
the wind, and W is the width of the large dune (Fig. [8]). 



The asymmetric dune shapes produced due to cohisions 
are different from the ones in Section [5l since the cohision 
not only imphes an asymmetric influx but also leads to a 
(dynamic) modification in the s hear stress fie l d dur ing the 
interaction between the dunes (jDuran et all l2QQ5l ). The 
hybrid, asymmetric dunes depicted in Fig. [9] are transient 
shapes. A sufficiently long time after the collision, dunes 
are well separated and adapt to constant influx and wind 
direction again. 
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Figure 8: Binary collision between barchans of different 
sizes. A small barchan (width ui) approaches a larger one 
(width W) from behind with lateral offset. The small and 
the large barchans are centered at y and F, respectively. 
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Figure 9: Snapshots of simulations of binary collisions for 
different values of A and volume ratio r ex {w/W)^ . Cal- 
culations with r = 0.06 are shown for (a) A = 0.2, (c) 
A = 0.6 and (e) A = 0.9; calculations with r = 0.3 are 
shown for (b) A = 0.2, (d) A = 0.6 and (f) A = 0.9. 



The simulations corroborate previous observations that 
dune asymmetry could result from limb merging as two 
barchans link laterally, ther eby caus i ng do wnwind exten- 
sion of the coalesced limbs (jBourkd , I2OIOI ). Such a situ- 
ation is found for sufficiently large offset (A > 0.9), as 
can be seen in Figs. [9^ and [9f. It is interesting to notice 
that the result of the lateral coalescence is the release of a 



small dune from the merged limbs, whereas the colliding 
barchans merge then into a single asymmetric dune. For 
other values of A, asymmetric barchan shapes can form 
during the separation of the barchans resulting from the 
collisional interaction. However, no limb extension occurs 
as the faster barchan upwind is approaching the larger one 
downwind. 

Asymmetric dune patterns very different from the ones 
in Fig.[9]can result from collisions between multiple dunes. 
Figure [TOk shows one snapshot of the simulation of a 
barchan dune field emerging from a flat sand surface sub- 
jected to a saturated influx. The simulations of the genesis 
of barchan fie l ds were discussed i n detail in p revious works 
(|Duran et al.L I2OIOI : iLuna et al.L I2OIII l2012f ) . As shown in 
these works, the small transverse bedforms emerging from 
the instabilities which develop on the flat hill upwind give 
place to small barchans after reaching the bedrock — due 
to the cross-wind instability inherent to transverse dunes 
(iReffet et al.UgOlOl : IParteli et all I2OIII : iMelo et al.L 12012 : 
iNiiva et al.L |2012[ ). The barchan dunes, which are sub- 
jected to a strong influx, grow in size during their migra- 
tion, thereby experiencing multiple collisions with their 
counterparts and producing complex, asymmetric dune 
shapes. Different patterns produced in the simulation in- 
clude limb extension of the upwind barchan, as we can see 
in Fig. [TOk . While the simulation of Fig. [TOk denotes a 
model for an immature dune field (which is at the early 
stages of its development), we have also performed sim- 
ulations using periodic boundary conditions in the wind 
direction such as to model the ste ady-state dynamics of a 
dune corridor (jDuran et al.L I2OIOI ). Such simulations also 
lead to asymmetric barchan shapes which cannot be ob- 
tained with binary collisions (cf. Fig. [TOb). 

The complex asymmetric shapes shown in Figs.fTOk and 
[TOb result from collisions between multiple barchans which 
have different sizes and offset values, and are further sub- 
jected to an asymmetric influx depending on the spatial 
distribution of dunes upwind. The systematic study of 
collisions between multiple dunes in a field is out of the 
scope of the present study. However, on the basis of our 
calculations, we can conclude that binary collisions can 
produce a rather limited range of asymmetric dune pat- 
terns. The only previously reported asymmetric pattern 
produced by a binary collision is the limb elongation due 
to dunes merging laterally (Fig. [9^,f). Other asymmetric 
patterns induced by collisions can result if the collisions 
involve multiple dunes during the evolution of a barchan 
dune field. As can be seen in Figs. [TOk and [TOb . colhsions 
between dunes lead to the formation of groups of barchans, 
in which upwind barchans releasing sand through their 
extended limbs act as source of sediment t o the dunes 
in the front, as previously noted by Tsoarl (|2QQ1[ ). In- 
deed, the complex patterns emerging from collisions be- 
tween barchans are transient dune shape s, as shown in 
Fig, p and in previous works (Schwammle and Herrmanm 
l2003nDuran et al.L[2QQ5[l2010[ ). Further modeling work is 
required in order to deepen our understanding of the colli- 
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Figure 10: Asymmetric barchans emerge during the evolution of a barchan dune field, (a) Simulation of the genesis 
of a dune field. A small, flat sand h ill subjected to saturated flux becomes a source of sand for a barchan dune fleld 



( Duran et al.luQloULuna et al.U2Q12[ ). Collision between multiple dunes lead to complex asymmetric dune morphologies 
which are not obtained in binary collisions (cf. Fig. [9]). (b) Complex asymmetric barchan shapes also occur in the 
simulation of a barchan corridor using periodic boundary conditions both in the direction perpendicular and parallel to 
the wind. Wind blows from left to right in both (a) and (b). Simulations were performed with r/rt ~ 3.9. 



sion dynamics of multiple dunes and of the resulting dune 
patterns. 

7. Conclusions 

In conclusion, we have investigated, by means of mor- 
phodynamic modeling of aeolian dunes, the respective 
roles of bimodal wind regimes, topography, influx asym- 
metry and dune collisions on the formation and evolution 
of asymmetric barchan dune shapes. The conclusions of 
our numerical simulations can be summarized as follows: 

• Bimodal wind regimes — the limb opposite to the sec- 
ondary wind elongates into the resultant transport di- 
rection if the divergence angle of the bimodal wind is 
obtuse and the ratio [Q2^w2]/[Qi^wi] — where Qi 
{Q2) and Twi (^^2) stand for the bulk sand flux and 
the duration of the primary (secondary) wind, respec- 
tively — exceeds 25%. These conditions are the same 
as the ones for oblique alignment of b edforms in dense 



sand beds (|Rubin and Hunterl . 119871 ). the asymmetric 
barchan being the corresponding morphology under 
low amount of available sand. 

• Topography — a barchan crossing a topographic rise 
can become asymmetric; the limb closest to the to- 
pographic rise elongates downwind. The migration 
velocity of a barchan that is on a tilted surface has a 



downhill component proportional to the tilting slope; 
the preferential limb extension is downhill rather than 
downwind. 

• Influx asymmetry — the side of the barchan subjected 
to the larger influx increases in volume, whereas the 
opposite limb elongates downwind. The typical asym- 
metric morphology is a barchan with an elongated 
arm (the one subjected to lower influx). The asym- 
metry is of transient nature, since the dune migrates 
laterally towards the region of higher influx. 

• Dune collisions — asymmetry can be triggered due 
to barchan collisions with lateral offset. Binary colli- 
sions with large offset can lead to extension of the limb 
resulting from the merging of the two limbs partici- 
pating at the collision. Collisions between multiple 
barchans in a fleld can trigger more complex asym- 
metric patterns. Some patterns produced in these 
simulations involve the elongation of one limb of the 
upwind dune, as reported from previous observations 
(lBourkel . [2QTQf ). 



In order to improve the quantitative assessment of asym- 
metric barchans, the model should be extended in order to 
accou nt for secondary flow effects at the dune lee ( Tsoari 
2QQ1I ). which may be relevant for the dynamics of limb 
elongation due to bimodal wind regimes or due to topo- 
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graphic ridges with sharp slopes (jBourkd , |2Q1Q[ ). Fur- 
thermore, reahstic simulations of bimodal wind regimes 
should account for multiple secondary wind trends in con- 
sistence with the complex sand roses of real dune fields 
( Frvberger and Deam Il979l ). The results of our calcula- 
tions are potentially useful for inferring local wind regimes 
or spatial variations either in topography or in availability 
of mobile sediments in planetary dune fields where asym- 
metric barchans occur. 
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